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In this paper the cross sections of the processes e^e~ — >■ H + "f,H = »?c, ??c, XcO, Xci, Xc2 are 
calculated. The calculation is carried out at the leading twist approximation of light cone formalism. 
Within this approach the leading logarithmic radiative and relativistic corrections to the amplitudes 
are resumed. For the processes e'^e~ — >■ 77c, JJc + 7 one loop radiative corrections are taken into 
account. It is also shown that one loop leading logarithmic radiative corrections calculated within 
light cone formalism for the processes under study coincide with that obtained by direct calculations 
of one loop diagrams within nonrelativistic QCD. 
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CN ; I. INTRODUCTION 

C , 

^ Theoretical approach to the description of hard exclusive processes, which is called light cone formalism (LC), is 

based on the factorization theorem Within this theorem the amplitude of hard exclusive process can be separated 

into two parts. The first part is partons production at very small distances, which can be treated within perturbative 
QCD. The second part is the hardronization of the partons at large distances. For hard exclusive processes it can be 
, parameterized by process independent distribution amplitudes (DA) . 
<~i ■ The production of the charmonium meson H in the process e+e^ H + j at B-factories, is the simplest example 
Qh' of hard exclusive process. One can assume that the energy at which B-factories operate is sufficiently large so that 
it is possible to apply LC. Another approach to the calculation of the cross section of this process is nonrelativistic 
(]J QCD (NRQCD) This approach is based on the assumption that relative velocity of quark-antiquark pair in 
, charmonia is small parameter in which the amplitude of charmonium production can be expanded. LC has two very 
important advantages in comparison to the NRQCD. The first advantage is that LC formalism can be applied for 
light or heavy mesons if DA of this meson is known. From, NRQCD perspective this means that LC resums whole 
J> ■ series of relativistic corrections to the amplitude under study. For NRQCD relativistic corrections are very important 
00 ' especially for the production of exited charmonia mesons. The second advantage is that within LC one can resum 
0^ leading logarithmic radiative corrections to the amplitude in all loops. The main disadvantage of LC is that within 
this formalism it is rather difficult to control power corrections to the amplitude. 
, Within NRQCD the process e+e^ H + "/ was considered in papers [4-6]. In paper Q this process was considered 
' at the leading order approximation in relative velocity and strong coupling constant. The authors of paper \S\ took into 
, account one loop radiative corrections. In addition to the radiative corrections the first order relativistic corrections 
■ to the process e+e^ — > ?7c + 7 were considered in paper @. 

The only process considered within LC is e+e" — >■ rjc + j 0j B- The main drawback of these papers is that the 
authors used very simple model of DA of the rjc meson, which doesn't take into account relativistic motion in this 
meson. Recently, the leading twist DAs of charmonia mesons have become the object of intensive study (ol-flTj. The 
study of these DAs allowed one to build some models for charmonia DAs, that can be used in the calculation of 
different exclusive processes. 

In this paper the leading twist processes e+e~ H + j will be considered. Using helicity selection rules p^| - |20j it 
is not difficult to show that at the leading twist accuracy the mesons with longitudinal polarization and the following 
quantum numbers H =^ Sq,^ Pi,^ P2,^ P3 can be produced. So, in this paper the following processes will be considered: 
e+e" H + 7, H = r]c, r]'^, XcO, Xci, Xc2- To calculate the cross sections of these processes the model of DAs proposed 
in papers [TTI - ll3l ITg} will be used. 

This paper is organized as follows. In the next section the amplitudes of the processes under consideration will be 
derived. Numerical results and the discussion of these results will be given in the last section of this paper. 
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FIG. 1: The diagrams that contribute to the processes e+e —^H + ^,H = t/c, ?7c, Xco, Xci, Xc2 at the leading order approxi- 
mation, in the strong couphng constant 



II. THE AMPLITUDE OF THE PROCESS e+e" ^ H + j. 

In this section the leading twist approximation for the amplitude of the processes e^e^ — > H -\- ^, H ~ 
Vc,v'cyXcOjXci,Xc2 will be derived. The diagrams that contribute to the processes at the leading order approxi- 
mation in the strong coupling constant are shown in Fig. [T]As was noted in the introduction, selection rules tell us 
that at the leading twist accuracy all produced mesons are longitudinally polarized. So, the polarization vectors of 
these mesons are proportional to the momentum of these mesons. 

To calculate the amplitudes and the cross sections of the processes involved one needs the expressions for the 
following matrix element of the electromagnetic current J^(0): {H{p)'y{k)\J^{0)\0). For the production of the 
longitudinally polarized iJcVctXcI mesons it can parameterized as follows 

(i/(p)7(fc)| J^(0)|0) = Fh e^.^a.^e^p'^k^, (1) 

where e'^ is the polarization vector of the final photon. It causes no difficulties to find the expression for the formfactor 
FH=ric.ri'a-Xci the leading twist approximation 

where Qc is the charge of c-quark, the definitions of the constants .fri^,ri'^,Xci ^^'^ the DAs <t'ri^,ri'^.,Xai (d m) can be found 
in the Appendix, ^ is the fraction of relative momentum of the whole meson carried by quark- antiquark pair, ^ is the 
characteristic scale of the process, s = {p + fc)^. 

The expression for the production amplitude of the longitudinally polarized XcOi Xc2 mesons can be written in the 
following form 

(iJ(p)7(fc)| Jp(0)|0) = Fh [{eq)k^ - {kq)e^ , (3) 

where q = k+p. The other designations are the same as were used in equation ([l} . The expression for the formfactor 
FH=xao,x,:2 tis-s the form 

The constants fxao,x<:2 ^-nd the DAs 4>xao,Xo2i^^ 1^) can be found in the Appendix. 
The cross section of the processes can be written in the following form 
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It should be noted that the matrix elements of the processes under study were taken at the leading order approximation 
in Mjj/s. The factor 1 — Mfj/ s in the cross section appeared due to the phase space of the final particles. 

The expression of the formfactor Fh depend on the DAs (/)}{{£,, fJ-) of the charmonia mesons. If infinitely narrow 
distribution amplitudes 0r)c,»;' ,xci (Cj m) = ^Xco.x^' /^) ^ ^'^'(0 substituted to formulas than 

NRQCD results for the amplitude will be reproduced [4]. If real distribution amplitudes (Ph{£,,IJ') are taken at the 
scale /i ^ rric, than formulas ([2]), (jlj will resum the relativistic corrections to the cross section up to 0{l/s'^) terms. 
To resum the relativistic and leading logarithmic radiative corrections simultaneously one must take the distribution 
amplitudes at the characteristic scale of the process /i ^ y/s. The calculation of the cross sections will be 

done at the scale fx = \/s/2. 

It is interesting to note that it is possible find the leading logarithmic radiative corrections at one loop level using 
formulas ©, dU without calculation of one loop diagrams. Applying the approach proposed in paper one gets the 
results 
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where rric is the pole mass of the c-quark, the definition of the NRQCD matrix elements {0)s, {0)p can found 
in paper Cf = 4/3. Note that in the above equations it was assumed that renormalization group evolution 
of the DAs begins at the scale /^o ~ and ends at the scale ^ ^ ^/s. At the scale /io ~ mc the DAs are 
'^»?c,r;e.xci (£j = '^xco,Xc2(C: m) = —^'(0- Leading order NRQCD results ([6]) coincide with the results obtained 

in paper Q . The one loop leading logarithmic radiative corrections for the F,,^ coincide with the result of paper Q ■ 
The one loop leading logarithmic radiative corrections for the Fjjcj ^Xcoi ^Xcu ^Xc2 coincide with the results of paper 



The result ([2]) for the production of the pseudoscalar mesons T]c,r]'^ can be improved since there exists expression 
for the one loop radiative correction to this amplitude [2l|, . This expression can be written as follows j2l| 
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In the above expressions it is assumed that the DAs are ^ even. 

It is instructive to take the limit of zero relative velocity of quark-antiquark pair and compare it to the NRQCD 
result [6] . At leading order approximation in relative velocity expression ([7]) becomes 
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The second term in equation ([S]) is due to renormalization group resummation of the leading logarithms in the DA. 
The last term is one loop radiative corrections to the hard part of the amplitude. The factorization scale separates 
long distance dynamic of the chamonium meson parameterized by DA from the small distance effects parameterized 
in the hard part of the amplitude. It is seen that /i dependence is canceled in the final answer, as it should be. 
The authors of paper [(| obtained the following NRQCD expression for equation ([S]) 
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It is seen that this expression is very similar to (|5]). Moreover one has one free parameter /io in expression 
can be used to adjust ([5]) to ([S]). However, expressions (|S]) and ^ seem to be a little bit different. 
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TABLE I: The cross sections of the processes e+e H + H = r/c, ??c5 XcO, Xci, Xc2. Second column contains the results 
obtained in this paper. In the third, fourth and fifth columns the results obtained in papers 0], ^] are shown. 



III. NUMERICAL RESULTS AND DISCUSSION. 



To obtain numerical results for the cross sections of the processes under study the following numerical parameters 
are needed. 

In this paper we are going to use the models of the charmonia DAs proposed in papers (iTI - fTsl [Tg} . For the strong 
coupling constant we use one-loop expression 
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where bo = 25/3 and Aqcd = 0.2 GeV. 

In the calculation the following values of the constants fu will be used [2S 



/xco(^J/*) 

fxcl 



0.373 ±0.064 GeV, 
0.261 ±0.077 GeV, 
0.093 ±0.017 GeV, 
0.272 ±0.048 GeV, 
0.131 ±0.023 GeV 
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The values of the constants /^^ , f^/ were calculated in paper [24] . The values of the constants of the P-wave charmonia 
mesons can be found in paper jlg| . It should be noted that the constants f^^^g , f-^^^ depend on the renormalization 
scale. As it is seen from formulas ([TU)) these constants are defined at the scale /j, — Mj/^,. The anomalous dimensions 
of these constants, which govern the evolution, can be found in paper (l6j . 

There are different sources of uncertainty to the results obtained in this paper. The most important uncertainties 
can be divided into the following groups: 

1. The uncertainty in the models of the distribution amplitudes 4>h{x, ^J-), which can be modeled by the variation of 
the parameters of these models . The calculation shows that this source of uncertainty is not greater than 10%. So, 
it is not very important and it will be ignored. 

2. The uncertainty due to radiative corrections. In the approach applied in this paper the leading logarithmic 
radiative corrections due to the evolution of the DAs and the strong coupling constant were resummed. For the 
processes e'^e~ — >■ r]c,r]'^ + 7 one loop radiative corrections were taken into account. So, for last two processes 
radiative corrections are not very important and they will be ignored. As to the other processes considered in this 
paper radiative corrections to the results can be estimated as as{s) ~ 20%. 

3. The uncertainty due to the power corrections. This uncertainty is determined by the next-to- leading order 
contribution in the 1/s expansion. One can estimate these corrections using the leading order NRQCD predictions 

as it was discussed in paper [5^. Thus, for the processes e'^e~ — >■ t?c, ??c, XcOj Xcij Xc2 + 7 the errors due to this 
source of uncertainty are ^ 3%, 6%, 50%, 37%, 60% correspondingly. 

4. The uncertainty in the values of constants ilO\) . The calculations show that, for the processes e+e^ — > 
Vcj v'cy XcO, Xcij Xc2 + 7 the errors due to this source of uncertainties are ^ 34%, 60%, 35%, 35%, 35% correspondingly. 

Adding all these uncertainties in quadrature one gets the total errors of the calculation. 

The results of the calculation are presented in Table Hill Second column contains the results obtained in this paper. 
In the third, fourth and fifth columns the results obtained in papers [!'], ^ are shown. It is seen that the results 
obtained in this paper are in reasonable agreement with the results obtained within NRQCD. 
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Appendix A: Distribution amplitudes. 



The leading twist distribution amplitudes needed in the calculation can be defined as follows: 
for the pseudoscalar mesons P = rjcrj'^: 

1 



for the Xco-meson: 



Xcoip) Qiiz)[z,-z]Q^ij{-z) 



= 



for the Xci-meson: 



Xcl(p,eA=o) Qaiz)[z,-'Z]Q^i^{-z) Oj = (p75)/3^ ^ g 



for the Xc2-meson: 

'Xc2(p,eA=o) Qiiz)[z,-z]Q^^{-z) o) = {p)pc 



The factor [z, — z], that makes the above matrix elements gauge invariant, is defined as 

[z,-z] = Pexp[ig dx'-' Af^{x)]. 

J —z 

It is not difficult to show that the functions (^,,^ ) and (j))^-^ (^) are ^-even. The normalization condition for these 
functions is 



1 



= 1. 



The functions 0xo(C) ^^'^ 't'x^iC) ^-re ^-odd and normalized according to 
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